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Abstract 

The (/-state clock model with the cosine potential has a single phase transition for q < 4 and 
two transitions for q > 5. It is shown by Monte Carlo simulations that the helicity modulus for 
the five-state clock model (q = 5) does not vanish at the high-temperature transition. This is in 
contrast to the clock models with q > 6 for which the helicity modulus vanishes. This means that 
the transition for the five-state clock model differs from the Kosterlitz-Thouless (KT) transition. 
It is also shown that this change in the transition is caused by an interplay between the number 
of angular directions and the interaction potential: by slightly modifying the interaction potential, 
the KT transition for q = 6 turns into the same non-KT transition. Likewise, the KT transition is 
recovered for q = 5 when the Villain potential is used. Comparisons with other clock-model results 
are made and discussed. 
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I. INTRODUCTION 



The research on two-dimensional (2D) continuous phase transition has a long history. 
One of the corner stones was Onsager's solution of the 2D Ising model [l| which led to 
the understanding of the connection between a broken symmetry and the universality of a 
continuous phase transition [1]. Another corner stone was the discovery of a 2D topological 
phase transition by Kosterlitz and Thouless (KT) 

a a. 

This topological transition does 
not involve any broken symmetry and has hence a completely different character. The q- 
state clock models have the interesting feature that both of these types of phase transitions 
are possible and to some extent interfere with each other. Partly because of this, the 
phase diagram for the g-state clock model is a longstanding question with as yet no clear 
consensus, both as regarding the number of phase transitions and the character of the phase 
transitions j^-ll]. The overall feature, for which there is full consensus, is that for q < 4, 



there is only one transition, whereas for q > q c there are two transitions. However jthe value 
of q c has been suggested to be both q c = 5 as in Refs. j^- 7] and q c = 6 as in Refs. [8, 9]. The 
most common view is that q r = 5 and that the two phase transitions for q > 5 are of the 

n q nn 

KT type [4K7|. If there is only one transition for q = 5 as suggested in Refs. |8|, |9J, it should 
be discontinuous More recently, it has been suggested that the upper transition for the 
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111 ] . Quite to the contrary, however, this 



121. 



six-state clock model is not the KT transition 
transition was later found to be of the KT type in Ref. 

In the present paper, we reinvestigate the transitions for q = 4, 5, and 6, using Monte 
Carlo (MC) simulations. It is found that the five-state clock model has two separate transi- 
tions. However, the upper transition is not a pure KT transition, since the helicity modulus 
does not vanish at the transition. It is also found that if the interaction potential is changed 
to the Villain potential, then the upper transition does become a pure KT transition. In 
a similar way, it is found that the six-state clock model has the KT transition and that it 
turns into a non-KT transition when the potential is slightly changed. If the potential is 
further changed the upper and lower phase transitions coalesce into one and become dis- 
continuous. Thus both the number of phase transitions and the character of the transitions 
for a clock model with q > 5 depend on the precise shape of the potential. This work is 
organized as follows: Sec. [TT] introduces the spin models and statistical quantities we are 
investigating. These models have somewhat different interaction potentials between spins 
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but have the same symmetry and quadratic leading term. We present our main numerical 
results in Sec. IIIII and summarize them in Sec. HVl 



II. CLOCK MODELS AND HELICITY MODULUS 

The g-state clock model is, like the Ising model, a model of interacting spins where the 
interaction energy between two neighboring spins is determined by the difference in spin 
angles. Hence we begin with a Hamiltonian derived from a spin-interaction potential U as 

H = Y t u{e i -9 j ) y (1) 

(ij) 

where 0$ is the ith spin angle, and the sum runs over all the nearest-neighbor pairs. We 
will here consider the case when the lattice is a 2D L x L square lattice. The interaction 
potential for the clock model is given by 

U((f>) = V((f>) = -Jcos(f>, (2) 

where J is the coupling constant. Each spin can only have q discrete directions, 9 = 
0^,..,^. This means that the g-state clock model has a discrete Z q symmetry. 
In case of the Ising model, q = 2, the possible spin directions are 6 = and tt, corresponding 
to the usual spin- up and spin-down directions, and the symmetry is Z 2 . In the limit q — > oo, 
the spin angles 9i becomes continuous and the Hamiltonian recovers a full U(l) symmetry. 
This limiting model is called the 2D XY model and is the prototype of a system exhibiting 
the KT transition. 

The critical properties of a continuous phase transition is to a large extent determined 
by the symmetry and the dimension of the model. Thus it is reasonable to assume that a 
small change in the interaction potential V(<p) will not matter as l ong as the Z q symmetry 



is preserved. Such a change is given by the Villain approximation [131 ]. 

U(4>) = Wmain(0) = -~ In | exp [-/3(<f> - 27m) 2 /2) I , (3) 

where is the inverse of temperature T in units such that the Boltzmann constant is unity. 
Figure [1] compares the clock- model potential V(<f>) with the Villain approximation. The 
point with the Villain approximation is that this approximation makes theoretical analysis 
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FIG. 1: (Color online) Interaction potentials in Eq. (Jl]) at various p values, compared with the 
Villain approximation, Eq. ([3]). 

more tractable j^l and it has often been tacitly assumed that the phase diagram and phase- 
transition properties obtained within the Villain approximation should be universal and also 
valid for the clock- model potential, Eq. 02]) [4j, |5|. 

In order to investigate the sensitivity to the precise functional form V(4>) in a systematic 



way, one can use the following parameterization 



2 7 

u(4>) = vm = ^ 



14 



15| 



cos 



2p 2 



(4) 



Note that V v= \ — J = — Jcos<p = V(4>). For p > 1 the potential has a smaller dip, and for 
p < 1 a larger dip, as is shown in Fig. [TJ Also note that for small 0, the leading term J0 2 /2 
is quadratic and is identical for V, Vyuiain; an d V p . 

In the present paper, we will focus on the helicity modulus which measures the resistance 
with respect to a uniform twist across the sample in one direction. For a twist of size A x 
across the x direction, this means that the Hamiltonian in the presence of the twist field 
is given by H(A X ) = Yluj) Vp{@i — % — A x /L x ), where L x = L is the system size in the x 
direction. The helicity modulus T measures the increase of the free energy F caused by the 
twist in the limit of small A x . This increment is given by ^^\^ x= qA 2 x /2 = TA^/2, where 
the leading dependence in A x is of the second order since a twist will always increase the 
free energy or leave it invariant. For an interaction potential U ((/)), the helicity modulus is 
given by 

T = (e)-L^(s 2 ), (5) 



<ii> : 



U' 1 



-9,) and s = L~ 2 u> 



9j) where the sum is over all 



with e = L~ 2 

links in the x direction. The derivatives are with respect to the argument so that U'(<f>) = 
dU/dcj) and U" = d 2 U/d<f) 2 . We will here also study a higher-order correlation function 
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FIG. 2: (Color online) 2D Ising model (q = 2). (a) Helicity modulus and (b) its temperature 
derivative, (c) The peak heights of S shows a logarithmic growth with the system size (see text). 
Error bars are shown in all of these plots. 

S = dT/dT = -l3 2 dT/d(3 written as 

5 = f { (eH) - (e) (H) + L 2 (s 2 ) - L 2 /3 [(s 2 H) - (s 2 ) (H)] } , (6) 
and check its size dependence H ~ L ah through MC simulations. 



III. MONTE CARLO SIMULATIONS 



We use the Wolff single-cluster algorithm 16| throughout this work. All the measurements 
presented here have been obtained by updating clusters O(10 6 ) up to O(10 7 ) times after 
equilibration. In Fig. EJ^a), we illustrate our calculations for the two-state clock model, 
which corresponds to the 2D Ising model. At the Ising transition, the helicity modulus 
makes a transition from a higher to a lower value. However, the helicity modulus remains 
positive and nonzero for all temperatures. This is in contrast to the KT transition for which 
the helicity modulus is zero in the high-temperature phase. The Ising criticality is instead 
reflected in the helicity modulus by the fact that its temperature derivative diverges at the 
critical temperature. This divergence is directly picked up by the correlation function S, 
as illustrated in Fig. MJo). From the peak heights in Fig. E£b), one can determine the size 
scaling E ~ L° h and the critical index ah- Figure [2|c) illustrates that = for the Ising 
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FIG. 3: (Color online) 2D XY model (q = oo). (a) Helicity modulus and (b) its temperature 
derivative, (c) The peak heights of E shows a divergence as L grows. Error bars are also shown. 

model since peak heights are consistent with a logarithmic growth. For the Ising model, this 
is a trivial result since all correlations involving s vanishes by s = 0. This mean that S is 
proportional to the specific heat for the Ising model and hence = a/u = 0. 

The other extreme for the clock models is the limiting case q = oo, which corresponds 
to the 2D XY model, the prototype model for the KT transition. Figure E](a) shows the 
helicity modulus in the 2D XY model for various sizes L. The characteristics for the helicity 
modulus in case of the KT transition is that it jumps from a finite value 2T c /7r to zero at 



T r in the limit L = oo 
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]. Figure [3]^a) is consistent with such a behavior. Figure EJ^b) 
illustrates that the peak heights of H diverges and the divergence is again shown to be 
consistent with a logarithmic growth [Fig. [3(c)] . The reason is, however, different since the 
correlation length £ for the KT transition diverges as ln£ oc (T — T c ) -1 / 2 jjj. At the same 
time, 5 is proportional to |T — T c | -1 / 2 [3], leading to the size scaling H ~ InL. Our MC 
simulations suggest that H for all q are consistent with at least a logarithmic size divergence 
[see Fig. IU(e)]. However, our MC simulations do not have enough precision at the largest 
sizes to rule out a stronger divergence in a power-law form. The crucial point in the present 
context is that the correlation function H displays a phase-transition singularity, manifested 
in a size divergence, for all the clock models including the Ising and the XY model. 

In Fig. HI we present the results for q = 3,4, 5, and 6. For each q, the helicity modulus T 
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FIG. 4: (Color online) Helicity modulus measured for clock models with (a) q = 3, (b) q = 4, (c) 
q = 5, and (d) q = 6, together with error bars, (e) The peak heights of S at the g values. 

and the higher-order correlation function S are given. Figure Hfe) shows that for all g, S has 
a critical divergence. This means that in all cases there exists a phase transition which can 
be associated with the helicity. An interesting point to note is that for q — 2, 3, 4, and 5, 
the helicity modulus itself remains finite for all temperatures. This means that in these four 
cases the transition is not of the KT type. The clock models with q = 2, 3, and 4 undergo 
well-known phase transitions: for q = 2 this is the Ising transition, for q = 3 the three-state 
Potts transition 2(3] and for q = 4 again an Ising-like transition 21|- These transitions 
go directly from the low-temperature to the high-temperature phase, which rules out the 
possibility for two consecutive transitions separated by a quasicritical phase characterized 
by a power-law decay of spin correlations. However, theoretical predictions strongly suggest 
that the five-state clock model does have two transitions 4-^]. This is also what we find 
from our simulations: Fig. Efa) determines the lower transition using the order parameter 
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FIG. 5: (Color online) (a) The order parameter in Eq. ([7]) around the lower transition temperature 
for the five-state clock model. Note that the crossing points between curves move slightly to the 
right as the size becomes larger, (b) Merging of Binder's cumulant around the higher transition 
point for the same model. 



introduced in Ref. [1 



= (cos(#)) , 



A8i 



me 



(7) 

itp 



where ip means the phase of the magnetization vector so that m = L ^\ e * 
This parameter distinguishes the true long-range order from the quasi-long-range order. On 
the other hand, Fig. E(b) detects the upper transition, at which the spin-correlations start 
to decay exponentially, using Binder's cumulant, 

|4\ 



u 



\m\ 



(8) 



2(|m| 2 ) 2 ' 

These measurements show that there are indeed two separated transitions. Thus the conclu- 
sion for the five-state clock model is that it does have two consecutive transitions. However, 
contrary to the theoretical expectations in Refs. 4j-|7| , it is not the KT transition since as is 
illustrated in Fig. HJc), the helicity modulus remains finite for all temperatures precisely as 
for q = 2,3, and 4. Nor is it a discontinuous transition since we do not observe any double 



peaks in the energy distribution 



15| 



In order to better understand the reason for this, we first note that the conclusions in 
Refs. were obtained using the Villain approximation. In Fig. El^a), we show the 

result for the helicity modulus for the five-state clock model using the Villain potential 
(compare Fig. [1]). In this case, the helicity modulus does indeed vanish. This means that, 
in accordance with Refs. {4- 7], the five-state clock model within the Villain approximation 
does have two consecutive transitions where the higher one is the KT transition. The crucial 
point made here is that this is not true for the real five-state clock model. The reason for 
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FIG. 6: (Color online) The helicity modulus for (a) the five-state Villain model, (b) the generalized 
six-state clock model with p = 2.0, and (c) the four-state Villain model. Error bars are shown but 
usually smaller than the symbol sizes. 

this change in the transition is clearly the interplay between the number of clock states, q, 
and the detailed shape of the potential. This is further illustrated in Fig. [6(b) , which shows 
the helicity modulus for the six-state clock model using the generalized potential in Eq. (jlj) 
with p = 2.0, which gives a slightly natter potential (Fig. [1]). For this potential, the six-state 
clock model has two consecutive transitions where the helicity modulus does not vanish at 
the upper transition, just as for the five-state clock model with the usual cosine potential. 
However, the six-state clock model has two consecutive transitions where the upper one 
is the KT transition both for the usual cosine potential [Fig. 11(d)] and for the Villain 
potential. Figure |6(c) shows the helicity modulus for the four-state clock model within the 
Villain approximation. However, in this case, the helicity modulus remains nonzero for all 
temperatures as for the usual cosine potential. 

IV. CONCLUSIONS 

From numerical simulations, it was verified that the five-state clock model has two transi- 
tions. It was also found that the helicity modulus T does not vanish at the upper transition. 
Yet it was found that the helicity modulus does vanish at the upper transition provided 
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the cosine potential is replaced by the Villain potential. In this latter case, the transition 
is perfectly consistent with the standard KT transition. Heuristically, this transition can 
be interpreted in terms of the unbinding of vortex-antivortex pairs: below the transition 
temperature, these pairs are bound together whereas above it the effective vortex-anivortex 
attractive interaction is softened such that free unbound vortices appear. It is the presence 



of such free vortices which is signaled by the vanishing of T [17H19| . It seems plausible that 
a description of the phase transition in terms of vortex-antivortex pairs will remain adequate 
even after the tiny change from the Villain potential to the usual clock-model potential. If so, 
this means that the rapid drop of T at the transition for the five-state clock model is again 
associated with a softening of the effective vortex-antivortex interaction, with the important 
difference that this softening is not enough to create free vortices. In this interpretation, it 
is the lack of free vortices which forces T to remain finite for all temperatures. Whether or 
not the softening will be enough to create free vortices then depends on an intricate balance 
between the number of clock states, q, and the precise form of the potential. This interpre- 
tation is consistent with the finding that, whereas the six-state clock model undergoes the 
KT transition for which the helicity modulus vanishes, it only takes a small modification of 
the potential using the systematic parameterization in Eq. (j3J) to change the transition such 
that the helicity modulus does not vanish. From this perspective, the topological transition 
of the five-state clock model, for which the helicity modulus does not vanish, appears to be 
a new and weaker cousin of the proper KT transition. 
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